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In this letter, we propose a general real-space method for the generation of nonparaxial accelerat¬ 
ing beams with arbitrary predefined convex trajectories. Our results lead to closed-form expressions 
for the required phase at the input plane. We present such closed-form results for a variety of 
caustic curves: besides circular, elliptic, and parabolic, we find for the first time general power-law 
and exponential trajectories. Furthermore, by changing the initial amplitude we can design different 
intensity profiles along the caustic. 


Since 2007, when exponentially truncated Airy beams 
were predicted [l| and observed [2[, the study of accel¬ 
erating waves has attracted a lot of attention (see for 
example the recent review article i). Such accelerating 
Airy waves follow parabolic caustic trajectories and were 
first predicted in the context of quantum mechanics Q- 
Furthermore accelerating beams with trajectories differ¬ 
ent than parabolic have been found The interest in 

accelerating waves is triggered by several associated ap¬ 
plications including filamentation (sf, particle manipula¬ 
tion micromachining ini , and Airy plasmons [H- 

In higher dimensions the radial caustic collapse of an 
accelerating beam generates a focusing effect abrupter as 
compared to the Lorenzian distribution of a lens focused 
Gaussian beam [T5|. Such abruptly autofocusing waves 
have found applications in generatin g a blation spots [T6| 
as well as in particle manipulation [^. Besides optics, 
curved and accelerated waves have been considered in 
different settings such as in electron beams 0 , in mat¬ 
ter waves [l8[ , and in antenna arrays 0. In addition to 
caustic beams, Bessel-like beams can be made to bend 
along arbitrary trajectories [^. [^. 

In the nonparaxial regime caustic beams offer the ad¬ 
vantage of bending at large angles. In this respect, 
accelerating beams following circular trajectories that 
can bend u p to an almost 90° angle have been demon¬ 
strated [^, . Such beams with elliptic and parabolic 

trajectories have also been found Nonparaxial 

accelerating beams following arbitrary trajectories have 
been observed in by numerically determining the 
Fourier space phase profile and in using plasmon- 
ics. Nonparaxial Bessel-like beams generated due to the 
presence of an extended focus, rather than a caust^ fol¬ 
lowing arbitrary trajectories have been proposed |^ . 

In this work we propose a general real-space method 
for the generation of nonparaxial caustic accelerating 
beams that follow arbitrary predefined convex trajecto¬ 
ries. More importantly, we present specific classes of tra¬ 
jectories where the phase is computed in closed form. 
Specifically, we analytically compute the initial phase 
profile required for nonparaxial accelerating beams fol¬ 
lowing not only circular, elliptic, and parabolic (which 
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are described in terms of known functions) but also for 
the first time for general power-law and exponential tra¬ 
jectories. By changing the initial amplitude profile we 
can generate accelerating beams with an intensity profile 
that is more uniform along the caustic or more intense in 
different locations of the trajectory. 

Let us consider the propagation of a monochromatic 
optical wave in a homogeneous isotropic medium with 
refractive index n. The divergence condition V • .E = 0 
is satisfied by choosing the electric field to be polarized 
along the ^-direction and to depend only on the trans¬ 
verse coordinate x and the propagation coordinate z as 
E = y'ip{x, Then the optical field amplitude sat¬ 

isfies the Helmholtz equation in two dimensions 

(V2 + fc2)^ = 0, (1) 

where k = nuj/c = 27r/A and = 5^ + 9^. Using 
the Rayleigh-Sommerfeld formulation the dynamics of an 
arbitrary initial condition '0o(O = is given by 
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where G{x,z]^) = —{i/A)H^\kr) is Green’s function of 
the Helmholtz equation, is a Hankel function, and 
r = \r\ = (x — + z‘^. Utilizing the asymptotic form 

HQ^\kr) rsj lA and applying the sta¬ 

tionary phase approximation we obtain the ray equation 
^ which after some calculations becomes 

x = C-k',{M0)z. ( 2 ) 

where 

(3) 


is the dispersion relation as obtained by selecting the first 
quadrant from k‘^ k‘1 = k‘^. 

Applying the implicit function theorem on Eq. (|2]) 
we derive the caustic trajectory of the beam as 1 = 
zk'^{(l)^{^))(l)^^{^).We would like to generate a nonparax¬ 
ial accelerating beam following a predefined convex tra¬ 
jectory of the form Xc = f{zc)- Noting that the rays 
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should be tangent to the trajectory, we obtain the follow¬ 
ing relation between the beam trajectory and the initial 
phase of the beam 

KiMO) = -fizeiO), (4) 

where the function Zc{^) is calculated from ^ = f{zc) — 
f'{zc)zc. Utilizing (|4]) we can compute from which 

the initial phase is derived by direct integration. 
This procedure can be implemented either analytically 
or numerically, depending on the choice of the function 
f{zc). In the following paragraphs we are going to apply 
this method to particular classes of trajectories in which 
the calculations are carried out analytically leading to 
closed-form expressions for the phase. 

As a first example, let us consider nonpar axial ac¬ 
celerating beams following circular trajectories mil. 
Rather than choosing a Bessel function as initial condi¬ 
tion, our method provides analytic expressions for the 
phase, allowing us to arbitrarily choose the amplitude 
function to amplify or suppress the intensity profile along 
specific parts of the trajectory. Selecting the first quad¬ 
rant of the circle 


4 + ( 5 ) 

results to Xc = f{zc) = ~ ? where R is the caustic 

radius. The relation between the transverse coordinate 
at the input plane and Zc is then given by 

zcio = wove - R^- ( 6 ) 

Combining Eqs. m, (© and m we find that = 

kx{^) = — R?/^. Finally, integrating the above 

equation, we obtain the following expression for the ini¬ 
tial phase 

m = kR[V{^lRf - 1 - sec -1 (^/i?)]. ( 7 ) 

The initial field amplitude is taken to be A{^) = 0 for 
0 < ^ < R and A{^) = for ^ > R, where 7 > 0 . 

Typical examples of accelerating beams with circular 
trajectories are shown in Fig. [TJ The simulation results 
are obtained by numerically solving Eq. O- Note that 
in all the figures the coordinates x, 2 ; are scaled accord¬ 
ing to the wavelength A. The effect of the choice of the 
initial amplitude function is presented in the top row; By 
changing 7 , we utilize the fact that rays emerging from 
smaller ^ contribute to earlier stages of the caustic tra¬ 
jectory (smaller values of the Zc) making the caustic more 
prominent at different locations. Specifically, in Fig.[T]^a) 
7 = 1/4 and the caustic is stronger for polar angles close 
to 0 = 7 r/ 2 . In Fig.[T]^c), 7 = 5 and the caustic is more in¬ 
tense closer to the initial plane. For intermediate values 
of 7 the caustic focuses in between as shown in Fig. [TJb) 
for 7 = 1. By engineering the initial envelope, the caus¬ 
tic remains quasi-invariant to an almost 90*^ arc as shown 
Fig. [IJd). The effect of the size on the beam trajectory is 
shown in Figs. [He)-(f) where the circle radius is reduced 
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FIG. 1. Amplitude profile of accelerating beams along circular 
caustic trajectories. In the upper panel R = 50A and (a) 
7 = 1/4, (b) 7 = 1, and (c) 7 = 5 . (d) R = 50A and 
MO = M0° + l/e + i/0R°, (e) R = 14A, 7 = 1/2 and 
(f) R = 7A, 7 = 1/2. The white/black dashed line represents 
the predefined trajectory of the beam. 


to R = 14A and R = 7X respectively. We note that, 
due to the generation of evanescent waves, the beam gets 
more distorted as the radius becomes smaller. 

Let us now consider accelerating beams following ellip¬ 
tic trajectories [234^ of the form 

Xc = fW = W - (8) 

Following the above described procedure we analytically 
computed the initial phase of the beam 

(j) = kR[E{sm-^ w(0k) - E{sin-^ w(-R)|r)], (9) 

where u{x) = x/{RW)^ r = 1 — and E{(p\m) is the 

elliptic integral of the second kind for a < 1 and 

(j) = kR[E{i sinh-^ u{R)\t) - E{i sinh"^ u{0\t)] (10) 

where u{x) = xj(RO— t), for a > 1. 

Typical examples of elliptic beam trajectories are 
shown in Fig.O Specifically, in the upper panel the caus¬ 
tics are constructed from rays launched from the major 
axis {a > R) of the ellipse, while in the lower panel the 
rays are launched from the minor axis {a < R) of the 
ellipse. The input wave envelope, A{^) has the same ex¬ 
pression as in the circular trajectory case. By changing 
the parameters of the amplitude function we generate in 
the three columns beams that are stronger in the ear¬ 
lier, the intermediate, and the later stages of the caustic, 
respectively. Such a flexibility might be useful for spe¬ 
cific applications, such as micromachining and particle 
manipulation. 

Let us now consider general power-law trajectories of 
the form 


Xc = f{zc) = azf 


( 11 ) 


where a > 0 and /3 > 1. Following the relevant calcula¬ 
tions we find that the input phase profile is given by 


m = 


aj3‘^kOaC)fi 

2^-1 


-1 


-2Ei{aMc\-e/^MV ( 12 ) 
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FIG. 2 . Amplitude profile of accelerating beams following 
elliptic caustic trajectories for R = 50A and a = 3R/2 (upper 
panel) and ol — Rj^ (lower panel), (a) 7 = 1/4 (b) 7 = 1, 
(c) 7 = 5, (d) 7 = 1/10, (e) 7 = 1/2, and (f) 7 = 5 . The 
white/black dashed line represents the pre-defined trajectory 
of the beam. 
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FIG. 3. Amplitude profile of accelerating beams following 
power-law trajectories for 7 = 1/2 and /3 = 3/2, (5 = 50A 
(upper panel) (a) d = 1, (b) d = 2, (c) d = 3; yd = 3 (middle 
panel) (d) d = 1, , (5 = 50A, (e) d = 1, (5 = 500A, (f) d = 3, 
8 — 50A ; yd = 5 (lower panel) (g) d = 1, (5 = 50A , (h) d = 1, 
8 — 500A, (i) d = 3, (5 = 50A. The white/black dashed line 
represents the analytical trajectory of the beam. 

where is a hypergeometric function, ^ = (yd —1)/(— ^) 

a = 1/2, 5 = 1 + l/[2(yd - 1)], and c = 2 + l/[2(yd - 1)]. 

The field amplitude in this case it taken to be A(^) = 
(d — (^)“^, where d > 0 for ^ < 0 and A(^) = 0 for 
^ > 0. If we scale the coordinates (xc^Zc) according to 
the wavelength, equation m becomes 

Xc/x = a{zc/X)^, ( 13 ) 

where d = is a dimensionless parameter. 

The upper panel of Fig. [3] depicts the amplitude profile 





of a power law beam trajectory with exponent yd = 3/2 
for different values of the scaling parameter d that af¬ 
fects the acceleration of the trajectory. As one can see, 
the beam amplitude is quite uniform along the caustic 
irrespectively of the value of d. In the middle panel of 
Fig. El beams with cubic trajectories (yd = 3) are shown. 
In Fig. [3je) we modify the amplitude function as com¬ 
pared to Fig.[3l^d) resulting to a more intense caustic for 
larger values of z. In Fig[3l^f) we change the scaling fac¬ 
tor d resulting to increased transverse acceleration. In 
the bottom panel of Fig. [3] accelerating trajectories with 
a larger exponent yd = 5 are presented. The presence 
of a “knee” in this case makes the intensity along the 
caustic trajectory less uniform and more intense around 
the “knee”. This happens due to the large amount of 
rays that contribute to the caustic around the “knee”. 
In Fig. [3l^h) we modify the initial amplitude profile as 
compared to Fig. ^g) resulting to a more uniform pro¬ 
file around the caustic after the “knee”. In Fig. ^i) we 
change d and thus the acceleration is increased. 

Finally, we apply the presented methodology for gen¬ 
erating beams that follow exponential trajectories having 
the form 


Xc = ( 14 ) 

We note that Xc{z = 0) = xoe~^^° = a and thus we 
should initially excite rays with ^ < (a. We select an ini¬ 
tial amplitude profile A{^) = (a + d — ^)“^ for ^ < a and 
A(^) = 0 for ^ > (a. Following the relevant calculations 
we compute the derivative of the initial phase 

d(j) _ (ayd/cexp (1 + FF (C)) 
di + (a^yd^ exp (2 (1 + FF (C))) ’ 

where FF is the Lambert FF function, C = —^/(ae), and 
e is the Euler’s number. Integrating the above equation, 
we find the initial phase profile: 

4>{0 = -9{0+9{a) ( 16 ) 

where g{C) = {k/j3) [tanh"^ /(C + (C)] , and 

fio = ^/T+we/w^■ 

In the upper panel of Fig. |4] accelerating exponential 
trajectories with different exponents are shown. As the 
exponent increases, the acceleration increases and the 
beam bends faster. In the lower panel of Fig. |l]we use 
different amplitude truncations leading to a caustic tra¬ 
jectory that is more intense in the beginning the inter¬ 
mediate or the later stages, respectively. 

Finally, let us point out that in all the figures one can 
see that the numerically obtained beam trajectories are in 
agreement with the expected geometrical curves (dashed 
white-black lines) [given by Eqs.(j5]), ([8]), (pTj) and ([T^ ]. 

In conclusion, in this letter we proposed a real-space 
method for generating nonparaxial accelerating caustics 
with arbitrary trajectories. We utilized the method for 
finding, in closed-form, the required initial phase for gen¬ 
erating accelerating beams following circular, elliptical, 
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FIG. 4. Amplitude profile of accelerating beams following 
exponential caustic trajectories for xq — 20A. In the upper 
panel zo = xo, S = xo and (a) [3 = 0.7/xo, (b) (3 — 3/2xo, (c) 
(3 — “3/xq. In the lower panel zq — 2xo, [3 = 3l2xo and (d) 
S = Xo, (e) (5 = 5xo, (c) S = lOa^o- In all cases 7 = 1/2. The 
white/black dashed line represents the analytical trajectory 
of the beam. 


and for the first time general power-law and exponen¬ 
tial trajectories. We would like to point out that similar 
methods can be applied to different nonparaxial settings, 
for example, for generating vector caustics or for de¬ 
signing different classes of nonparaxial trajectories [3l|. 
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